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The charge ordering phenomena in quasi two-dimensional 1/4-filled organic compounds 
{KT)2X (ET=BEDT-TTF) are investigated theoretically for the 9 and a-type structures, based 
on the Hartree approximation for the extended Hubbard models with both on-site and intersite 
Coulomb interactions. It is found that charge ordered states of stripe- type are stabilized for the 
relevant values of Coulomb energies, while the spatial pattern of the stripes sensitively depends 
on the anisotropy of the models. By comparing the results of calculations with the experimen- 
tal facts, where the effects of quantum fluctuation is incorporated by mapping the stripe-type 
charge ordered states to the S = 1/2 Heisenberg Hamiltonians, the actual charge patterns in the 
insulating phases of 0-(ET)2MM'(SCN)4 and a-(ET)2l3 are deduced. Furthermore, to obtain 
a unified view among the 9, a and k-(ET)2^ families, the stability of the charge ordered state 
in competition with the dimeric antiferromagnetic state viewed as the Mott insulating state, 
which is typically realized in K-type compounds, and with the paramagnetic metallic state, is 
also pursued by extracting essential parameters. 

KEYWORDS: charge ordering, organic conductors, BEDT-TTF, extended Hubbard model, Hartree approxima- 
tion. 



§1. Introduction 

The family of quasi two-dimensional (2D) organic conductors (ET)2X (ET=BEDT-TTF) is 
known to exhibit a variety of interesting electronic properties. & Their structure consists of alter- 
nating layers of anionic X~ with the closed shell, and cationic ET^/^+ whose vr-band is 3/4- filled 
(1/4-filled in terms of holes) if all the ET molecules are equivalent. The variety in physical prop- 
erties reflects that of spatial arrangements of ET molecules in the layer, which are classified by 
greek characters such as k, a, 0, /?, etc., together with the rather strong mutual Coulomb inter- 
action among 7r-electrons.& Theoretical studies of Kino and FukuyamaSi'i'& (KF) developed a 
systematic way to understand the diversity in their ground state properties by taking into account 
the explicit anisotropy of the transfer integrals between ET molecules for each compound, and by 
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mean-field (MF) calculations treating the on-site Coulomb interaction U within the Hartree-Fock 
approximation. 

By studying the structure of k-(ET)2X where the intradimer transfer integral is twice larger 
than the others,& and by varying the degree of the dimerization, KF have investigated the stability 
of the insulating state with antiferromagnetic (AF) ordering between dimers, which is called the 
dimeric AF state.& This state is stable when the dimerization in the FT layer is strong, since the 
electron correlation results in a Mott insulating state due to the 1/2-filled 7r-band,i'll' as has been 
first argued by Tokura© 

Another interesting conclusion of KF is that the a-type compounds show tendency toward an 
insulating state with charge ordering (CO) (charge transfer in their notation) S Experimentally, 
Q:-(ET)2l3 shows a metal-insulator (MI) transition at Tmi=135 kS-* at ambient pressure, and by 
applying pressure this Tmi is supressed.lli* The magnetic susceptibility shows a sharp drop below 
Tmi indicating a non-magnetic ground state with a spin gap.El) A bond length analysisB* as well 
as an IR absorbtion experimentlll) in this compound suggest charge disproportionation among FT 
molecules, though the suggested spatial patterns of charge disproportionation are different from 
each other. The Hartree-Fock calculations by KF for the structure of Q-(ET)2l3 show CO between 
the stacks I and II, which are the two crystallographically inequivalent columns as shown in Fig. 
1(a). The actual pattern of CO as well as the nature of the phase transition in this compound is 
still an open problem. Recently some suggestions have been made that the ground state of another 
class of a-phase, a-(ET)2MHg(SCN)4 for M=K, Rb, and Tl, which are isostructural to a-(ET)2l3 
(see Fig. 1(a)), may also be a charge disproportionated statelll^ rather than the so far believed 
spin-density- wave (SDW) state due to the nesting of Fermi surface.!' 

Recently such CO has been directly observed in NMR measurements by Miyagawa et a and 
by Chiba et aR) in a member of another polytype 9-(ET)2X . A series of these compounds, 9- 
(ET)2MM'(SCN)4 (M=Rb, Tl, Cs, M'=Co, Zn), has been synthesized by H. Mori et al^ who 
have proposed that the properties of 0-phase can be summarized in one phase diagram, where 
the horizontal axis is the dihedral angle, between the molecular planes, which is different 

for each compound. It has been argued by Takahashi and Nakumara© that a phase diagram, a 



slightly modified version of ref. 21, can be drawn as shown in Fig. 2, since it has been realized 
by experiments that among the above series, located in the region of 100° ^ </> ^ 140°, there exist 
compounds with two different types of structures, the usually called 0-type one, as shown in Fig. 
1(b), and another one with slight dimerization along the c-axis, as shown in Fig. 1(c), which we call 
here ^^-type to make distinction. The electronic properties of these two phases are distinct to each 
other. The compounds in the 0-phase which exist in the left side of the phase diagram (0 < 110°) 
are metallic at room temperature, and show a minimum in their resistivities at Tp = 50 — 100 K 
and by applying pressure Tp increases, i) except in 0-(ET)2l3 exhibiting superconducting behavior 
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Fig. 1. Schematic representation of the structures in the donor plane for {ET)2X compounds; a-type (a), 6-type 
(b) and 8d-type (c) (see the text for the notation of 9^)- ET molecules are drawn as ellipses, and the intermolecular 
transfer integrals are given by t = Es with E = —10 eV, where the overlap integrals s are for a-(ET)2l3 
where the indices are modified from those in ref. |l8| as al — a3 cl — c3 and 61 — 64 ^ pi — p4, (those 
corresponding to a-(ET)2NH4Hg(SCN)4 are also shown in the parentheses of (a)0), 6i-(ET)2CsZn(SCN)4B (b) 
and 6i-(ET)2RbCo(SCN)40'(c). 



below 3.6 K. Their magnetic susceptibilities show Curie-like behavior with no sign of magnetic 
order down to low temperatures.i3@) On the other hand, in the ^(^-phase located in the right side 
of the phase diagram {(p ^ 110°), the resistivity shows an insulating behavior {dp/dT < 0)0' 0) 
below room temperature and its magnetic susceptibility follows that of the ID Heisenberg model 
(J ~ 160 K for 0rf-(ET)2RbZn(SCN)4) but exhibits a spin gap behavior at low temperature in some 
compoundsBili) This spin gap behavior is also confirmed by ^^C-NMR measurements. Ej'EM) The 
CO was recently found in the ^f^-phase of (ET)2RbZn(SCN)4lll0 which is located in the critical 
region between the two phases.0 At room temperature it has ^-type structure, and upon cooling 
the sample slowly, the resistivity suddenly increases at T^tr = 195 K where it is believed that a 
structural phase transition takes place from 9 to ^^-type structure, while this structural phase 
transition is absent if the sample is cooled rapidly so that the 0-type structure is maintained down 
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CsZn RbZn Cu2(CN)[N(CN)2]2 
X= I3 CsCo RbCo TlCo 
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Fig. 2. Schematicphase diagram of 6-(ET)2X deduced from experiments, on the plane of dihedral angle {(j>) and 
temperature (T).0'0 The effect of pressure is indicated by the arrow, and SC denote the superconducting phase. 
The colored area show the S^-phase, while the other is the 0-phase. The locations of ambient pressure are shown 
for examples of anions X, where X = MM' (SCN)4 (M=Rb, Tl, Cs, Af'=Co, Zn) is abbreviated as MM'. 



to low temperatures. In the slowly cooled case, the ^'^C-NMR measurements show the existence of 
the charge disproportionation among ET molecules in the ^^-phase below Tstr, i.e. the occurence 
of a CO state.00) 

As for theoretical studies on the CO phenomena in 1/4-filled organic conductors, most are limited 
to those in one-dimensional (ID) cases. The effects of intersite Coulomb interactions Vij in such 
quasi ID systems have been studied from early dayli) and it is known that Vij makes the system 
toward the 4/cf CDW, i.e. CO, which is actually observed experimentally in (DI-DCNQI)2Ag. 
In the case of (TMTTF)2X the interplay between this Vij and the on-site Coulomb energy U 
together with the intrinsic dimerization give rise to a intermediated state between the dimeric AF 
state and the CO state,!!!) and in the cases of X =Br and SCN experiments suggest the existence 
of the expected charge disproportionation.!!® As regards the 2D case, the importance of this Vij 
has been clarified recently in a 1/4-filled transition metal oxide NaV205, where a similar interplay 
between U, Vij and the dimerization is present, and its ground state is proposed to be a peculiar 
type of CO state due to this Vij.© As for the organic systems, the calculations without Vij for 
6'rf-type structure predict the ground state to be Mott insulating state, 

i) 

inconsistent with the 

existence of CO disclosed by the experiments mentioned above. Thus, studies including Vij are 
needed to understand CO phenomena in the above 2D organic conductors, although CO itself has 
already been suggested in a-type structures from the Hartree-Fock calculations by KF 
The actual values of this Vij in ET compounds are expected to be relatively large due to the 
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fact that the wave function of the MO in molecular conductors is rather extended, so U is expected 
to be supressed so that to be comparable to neighboring Vi,j. Furthermore, the screening of Vij 
cannot be expected in ET compounds for the HOMO vr-band which is seperated from the other 
bands, as in the transition metals where the s-band screens the long range Coulomb interactions 
between the d-electrons. Actually several calculations for organic compounds estimate the nearest- 
neighbor Coulomb energy V to be 20 - 50% of ^7.00) In this paper we consider the value of U 
on ET molecule to be approximately 0.7 eV, which is derived from the analysis of NMR data of 
K-(ET)2Cu[N(CN)2]Br based on random phase approximation, and the values of the neighboring 
Vij to be 0.15 — 0.35 eV, adopting the ratio Vij/U = 0.2 — 0.5 from the above calculations. 

In this paper, we investigate the effect of both the on-site and the intersite Coulomb interactions 
on the electronic states in (ET)2X and our goal is not only to understand the electronic properties 
of the 9, 9(1 and a-phases but also to provide a unified view on the variety of ground states of 
this family. In §|2| the interrelationship between the 9, 9^, a and K-type structures based on a 
simplified model is explained together with the formulation of the MF calculations based on the 
work of KF,^! and also an explanation of the CO states of stripe-type, which is an analogy to the 
phenomena first found in Nickel oxides J^l* is given. The results of calculations for each type of 
structure, i.e. 9, 9^ and a-type, are shown in together with the results for the simplified model. 
The relevance to the actual compounds is discussed in §Q by comparing them with the experiments 
and a unified view on the electronic properties of (ET)2-'^ is given. §^ is devoted to the summary 
and the conclusion. In the following sections we call the a-type structures with transfer integrals 
corresponding to a-(ET)2l3 and a-(ET)2MIIg(SCN)4, as als-type and aMHg-type, respectively, 
for abbreviation. 

§2. Model and Formulation of Mean Field Calculations 

2.1 Interrelationship among different polytypes of (ET)2X 

It is useful to compare the structural features and the band structures of different polytypes, which 
has been carried out by KFi'S who pursued the interrelationship between k, als and aMHg-type 
structures on the plane of two key ingredients, i.e. the degree of dimerization and that of band 
overlap. Here we review their framework and also find out where the 9 and 0(i-type structures are 
located, by referring a simplified model described in Fig. 3. This model has been introduced by KF 
in order to unify their results for the a-type and «;-type structures, yet they have not performed 
actual calculation on this simplified model. 

In Fig. 4, the band structures of 0-type (a), 9^ and K-type (b), and als and aMHg-typeii' (c) 
are shown, respectively. It can be seen in Fig. 4(a) that the energy band of 0-type compounds 
is simply 3/4-filled, while in Fig. 4(b) for 9fi and K-types the upper two bands and the lower two 
bands are split due to the dimerization. Especially, the band of K-phase has a clear dimerization 
gap while the splitting is not appreciable in the ^^-phase, since the relative value of the interdimer 
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transfer integral compared to other transfer integrals, i.e. the degree of dimerization, is quite larger 
in the «;-type than in the Od-type; the interdimer transfer integral in the K-type structure is 0.257 
eV (tfti in Fig. |^ which is twice larger than the other transfer integrals while the one for the Od- 
type structure is 0.144 eV {tps in Fig. 1(c)), slightly larger than the other transfer integrals in the 
transeverse direction, tpi, tp2 and tp4. The difference among the two a- type compounds has been 
clarified in the work of KF,i'i) who indicated that this is the degree of the band overlap resulting 
from the difference in the absolute value of the transfer integral tp4 as shown in Fig. 1(a). The 
relatively small value of |tp4| in the als-type structure results in the small band overlap, while the 
energy bands for the aMHg-type structure has larger overlap due to the larger value of |tp4|, as 
seen in Fig. 4(c). 

It turns out that the simplified model in Fig. 3 is also an effective model for these k and Od- 
type structures by varying tpi, and for a- type structures by varying tp4, as can be seen in the 
following. The model in Fig. 3(a) is the 9-type structure which has the inversion as well as the 
spiral symmetry. If the absolute value of tpi in Fig. 3(b) is increased, the inversion symmetry is 
lost, and only the spiral symmetry is retained, which corresponds to the situation of k and ^^-type 
structures. If tp4 in Fig. 3(b) is varied, then only the inversion symmetry is retained, which is 
the situation of a-type structures. In Fig. 5, the band structures for the model in Fig. 3(b) are 
shown for different parameters, which are ta = 0.01 eV, tbi = h = 0.1 eV and thi = 0.1,0.2,0.3 
eV (a), and ta = 0.01 eV, thi = 4 = 0.1 eV and t^ = 0.1,0.05,0.0 eV (b). The general features 
seen in Figs. 4 (b) and (c) are well reproduced in Figs. 5 (a) and (b), respectively. Thus the 
calculations on this simplified model will provide us the idea of a unified view of the (ET)2X 
family. The interrelationship between the phases can be summarized in Fig. 6, which is similar to 
the corresponding figure in ref. ^ but the locations of the and 6'rf-phases beeing added. 




Fig. 3. A simplified model for (ET)2X introduced by KFo* The circles and the arrows represent the ET molecules 
and the intermolecular transfer integrals, respectively. 
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Fig. 4. The tight binding band structures of 9- type (a), ^^-type and K-type (b), ala-type and QMHg-typJi)(c) 
structures, where the transfer integrals are taken from those in Fig. 1. The fermi energy for each case is also 
shown. 



2.2 Formulation 

2D layer of ET molecules is considered by taking into account the transfer integrals betweeen them 
and the mutual Coulomb interaction among electrons. In this work the treatment of KF is basicly 
followed,!' but here the intersite Coulomb interactions are also included. Thus our Hamiltonian is 
an extended Hubbard one, which is as follows; 



H= Y (tijataja + h.c. 



<t,j> <^ 
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Fig. 5. The band structures for the simpUfied model shown in Fig. 3 (b). tc and tp are fixed at 0.01 eV and 0.1 
eV, respectively. The band structures for the case of tp4 = tp varying tpi = 0.1, 0.15, 0.20 eV (a) and for the case 
of tpi = tp and tp4 = 0.1, 0.05, 0.0 eV (b) are shown together with the fermi energy for each case. 




Fig. 6. Interrelationship between different types of structures for (ET)2X, schematically shown in the plane of tpi 
and tp4, which are the transfer integrals in the model structure described in the inset (identical to the one in Fig. 
3(b) ), which controls the degree of dimerization and band overlap, respectively. 
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i <i,j> 



(1) 



where < i,j > denotes the neighboring site pair, a is the spin index which takes | and J,, rii^ and 
(fljcr) denote the number operator and the creation (annihilation) operator for the electron of 
spin a at the ith site, respectively, and nj = nj| + 

The Coulomb interactions U, Vij are treated within the MF approximation in a manner similar 
to that in refs. |-|,|l|,|4; 



mn 



1,0 J 



(nil) mi + mi imi) - {nn) {mi) 
{m) nj + m. {nj) - im) {nj) , 



(2) 



which corresponds to the Hartree approximation. The MF Hamiltonian in /c-space is given by 



H 



MF 



E 

kcr 



0-2k(T 



ho + U 



( 



V 



n2a 







m 



+ hi 



J 



I 



a2kcr 



\ 



(3) 



\ ^mka ) 



where /iq and hy denote the matrices for the kinetic energy term and the intersite Coulomb 
interaction term, respectively. The actual form of these matrices depends on each structure and 
unit cell size whose examples are shown in § p.3| . 

The crysatallographic unit cell in the layer for the 0-type structure is approximated as rectangle 
with a = 2c and those for the Qd and a-type structures as square. 

Self-consistent solutions are searched for at T = with the average electron density being fixed 
at 1.5 per ET site by considering several unit cell sizes and their energies are compared so that the 
ground state can be decided. The solutions where the net magnetic moment of the whole system 
is finite, i.e. ferro/ferrimagnetic solutions are neglected, as in the work of KF, since the stability of 
these states, if it occurs, can be thought to be the artificat of the approximation. The total energy 
E is calculated as 

1 



A^cell 



E ^Ikanp {eiko 



Ika 



E im) (rij) , 

<id> 



(4) 



where N^eii is the total number of the cells, eika is the Ith eigenvalue of the m x m Hamiltonian 
matrix in eq. 3 at each k and np (e) is the Fermi distribution function. 

The intersite Coulomb interactions we consider are taken along the bonds with transfer integrals 
in Fig. 1, i.e. those between the neighboring ETs. For simplification only two kinds of values are 
considered: the intersite Coulomb energy along the stacking direction Vc which is the one between 
parallel molecules, and along the bonds in the transverse directions Vp which is the one between 
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molecules with dihedral angle (p. In the actual compounds, the values of intermolecular Coulomb 
energy should be different for the bonds with different transfer integral since the configurations 
between pair of molecules are different. So the above model with only Vc and Vp should be appro- 
priate for the 0-type structure with only two kinds of values of transfer integrals tc and tp. However 
we approximate in the calculations also for the 6d and a-type structures that only two kinds of the 
intersite Coulomb energy Vc and Vp exist since the deviation in the actual 6d and a-(ET)2X, which 
both can be considered as structurally similar phases to the 0-phase,© is expected to be small. 

2.3 Stripe-type charge ordered state 

In the following sections the solutions with charge disproportionation of stripe-type are found 
to be relevant in some compounds, where the sites with more amount of charge are arranged in 
rows alternately. Three kinds of charge patterns, namely arrangements of holes, can be naturally 
considered as shown in Fig. 7, which we call the vertical stripes (a), the horizontal stripes (b) and 
the diagonal stripes (c), respectively. In Fig. 7, the schematic representation of the solutions with 
stripe-type CO which are found to be stabilized are shown. In Fig. 8(a) and (b) two examples 
for the CO state with vertical stripes for the ^^-type structure are displayed where the unit cells 
consist of 4 and 8 ETs (m = 4 and m = 8 in eq. 3), respectively, where the matrix elements of /iq 
and hy for the former case of Fig. 8(a) are written as 



[^0]i3 

[^o]23 
[^o]24 



2tp2e-'''y'^cos{k^/2), 
2tpie-'''yl^cos{k,/2), 







([N34) 



tcie' 




[V]ii = 2Vcn2 + 2Vp{n3 + n^) 
[hv]22 = 2^471.1 + 2Vp{ns + 71.4) 
[/iy]33 = 214714 + 2Vp{ni + 712) 
[hvU = '^^cfis + 2Vp{ni + 712) 
[hvh = (i^j). 
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(a) vertical stripes 




(b) diagonal stripes 
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Fig- 7. Schematic representation of the spin moments and hole densities for the stripe-type CO solutions found to be 
stabilized. The diameters of the gray circles and the arrows represent the hole densities 2— < m > and the spin moments, 
respectively. 




Fig. 8. The magnetic cell for two examples of vertical stripe-type CO solution for ^d-type structure. Thick lines 
represent the bonds with the largest transfer integral tps. The dotted areas, the larger gray circles and the arrows 
denote the unit cells, the sites with excess holes and the magnetic moments, respectively. 
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§3. Results of Calculations 

In self-consistent MF calculations, the 2D integrations are approximated by a sum of the mesh, 
where the analytical tetrahedron method is employed!'© for 30x30 A;-points per Broullian zone 



for the cell of four sites. In §3.1 - B^, the value of U is fixed at U = 0.7 eV. The dependences on 



U are considered in S3.4 



3.1 e-type 

The results for 0-type structure reveal that stripe-type CO states have the lowest energy for 
relevant values of Coulomb energies, though the charge pattern depends on the parameters of 
the model, which are the transfer integrals tc and tp and the intersite Coulomb energies Vc and 
Vp. First we consider the case of tc = 0.01 eV and tp = —0.1 eV, the typical values for the 
series of 6'-(ET)2MM'(SCN)4,E3) and the case of the isotropic intersite Coulomb interactions, i.e. 
Vc = Vp = V. In Fig. 9, we show the ^-dependence of the calculated relative energies for the 
solutions which are stable for < V/U < 1 compared to the energy of the diagonal stripe solution 
schematically shown in an inset. The solutions which are not shown in the figure have higher 
energies. It can be seen that the vertical stripe solution with AF ordering both along and between 
these stripes (see inset of Fig. 9), has the lowest energy for V/U < (y/U)cTi = 0.49, whereas for 
{V/U)cri < V/U < {V/U)ct2 = 0.51, the state with diagonal stripes does. For V/U > {V/U)cv2 a 
nonmagnetic insulating state is stable, which can be seen as a bipolaronic state with every fourth 
site being occupied by almost two holes, as schematically shown in an inset. 

The value of {V/U)cri slightly varies when the values of transfer energies, tc and tp, are var- 
ied, though the qualitative behaviors are similar for tc/\tp\ < 0.2, which is appropriate for 9- 
(ET)2MM'(SCN)4© In Fig. 10 the phase diagram for the O-type structure is shown on the 
plane of V/U and tc, with fixed values of U = 0.7 eV and tp = —0.1 eV, excluding the bipolaronic 
state to see the competition between different stripe-type CO states. It can be seen that, concern- 
ing the charge pattern, the vertical stripes are preferred in the region of small V while the diagonal 
stripes are in the region of large V as an overall feature, though the spin configurations depend on 
the value of tc- 

Next we vary the value of Vp/Vc yet fixing U = 0.7 eV, tc = 0.01 eV and tp = —0.1 eV. In 
Fig. 11, the T^/T4-dependence of the relative energy of the vertical stripe solution compared to the 
diagonal stripe solution for the case of Vc fixed at 0.25 eV is shown. Note that the other solutions 
have higher energies than these two solutions. As Vp/Vc is decreased, there exists a phase transition 
at {Vp/Vc)cr = 0.95 from the CO state with vertical stripes to the one with diagonal stripes. It 
is noteworthy that a very slight deviation from Vp/Vc = !> where the state with vertical stripes 
is stable, makes the energy of the diagonal stripe solution lower, so the boundary betweeen the 
vertical and diagonal stripes in Fig. 10, which is around V/U = 0.5 for the isotropic case, is easily 
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Fig. 9. V lU dependence of relative energy compared to the diagonal stripe solution which is schematically shown 
in an inset, for \J = 0.7 eV, tc = 0.01 eV and tp = —0.1 eV. The schematic representations of the charge pattern 
as well as the magnetic structure in these solutions are shown in insets. The enlargement of the critical region is 
shown in (b). 



varied when Vp/Vc is varied. 

3.2 Od-type 

The calculations for the Od-^JV^ structure show that several different types of charge patterns 
compete with each other when the values of intersite Coulomb energies are varied. In this subsection 
the transfer integrals are fixed to those in Fig. 1(c). In contrast to the ^-type structure above, there 
exist three different charge patterns for horizontal stripe solutions which have different energies: 
the solutions with excess holes on sites 1 and 3 (stripes along tp3 and tpi), sites 1 and 4 (stripes 
along tp4), and sites 2 and 3 (stripes along tp2) as shown in Fig. 12(a)~(c), respectively. Among 
these states, the stripes along tp-j, and tpi have lower energies in the region of large V {V/U > 0.5), 
whereas for small V {V/U < 0.5) the stripes along tp4 does. In Fig. 13, the ground state solutions 
are shown for the case of isotropic intersite Coulomb interactions, Vc = Vp = V, as V/U is varied 
with U fixed at 0.7 eV. It can be seen there that when V is increased the charge pattern of 
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Fig. 10. The phase diagram for the d-type structure on the plane of V/U and tc, for U = 0.7 eV and tp = —0.1 eV, 
excluding the bipolaronic state. 
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Fig. 11. Vp/Vc dependence of relative energy compared to the diagonal stripe solution which is schematically shown 
in the right inset, on the ^-type structure for the fixed value of Vc = 0.25 eV and U = 0.7 eV. 



the solution having the lowest energy varies from the vertical stripes to the bipolaronic state at 
V/U = {V/U)cr = 0.50. In Fig. 14, the magnetic moment and the charge density on each site 
for the vertical stripe solution are shown as a function of V/U. As V is decreased, this stripe-type 
CO state continuously crosses over to the dimeric AF state, with almost same charge density on 
each site and AF ordering between dimers, i.e. pairs of molecules connected by the largest transfer 
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Fig. 12. Schematic representations of the charge pattern the three different horizontal stripes in 9d-type structure: 
stripes along tp3 and tpi (a), those along tp4 (b) and those along tp2 (c). The thick dotted, double dotted, thick 
solid, double solid bonds represent the bonds with the transfer integral tpi, tp2, tps, tpi, respectively and the thin 
solid bonds are for tci-2- 



energy tpz, which is known to be stable for V = 0.^3^ 

When Vp/Vc is varied, the solution which gives the lowest energy easily changes as in the case of 
^-phase. In Fig. 15, we show the variation of ground state solutions for the fixed value of 14=0.25 
eV. For Vp/Vc < (^/V'c)cr = 0.88, the state with diagonal stripes is stable, whereas there exists 
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Fig. 13. The sequences of the solution which gives the lowest energy for 0d-type structure, as a function of V/U 
{U — 0.7 eV). The thick bonds represent the bonds with the largest transfer energy tp^, as in Fig. 8. 
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Fig. 14. V/U dependence of the charge densities (a) and magnetic moments (b) on each site for the state with 
vertical stripes which gives the lowest energy in the region of small V for the 0d-type structure. Inset shows the 
indices for sites and schematically the charge pattern and the magnetic structure in this vertical stripe solution, 
where the thick bonds describe the bonds with largest transfer energy tpz, the sites with large (small) arrows and 
those with large (small) diameters represent the sites with large (small) hole density and those with large (small) 
magnetic moment while the direction of the arrows shows the direction of spins. 



another solution with almost same energy: the horizontal stripes along fp4, as shown in Fig. 15. 
For Vp/Vc > {yp/Vc)cT: the vertical stripe solution has the lowest energy, as in the case of 0-phase. 

3.3 al^-type and aM Hg-type 

The calculations for the a-type structures show that the diagonal stripe solutions have rather 
high energies compared to the other stripc-typc CO states, which is in contrast to the above cases 
for the 9 and Od-typc structures. In this subsection the transfer integrals arc fixed to those in 
Fig. 1(a). Due to the lower symmetry than the 0-typc structure, as in the fl^^-phasc, concerning 
the vertical and the horizontal stripes each has two solution with different charge patterns having 
different energies: as for the vertical stripes those along stack I and those along stack II, and for 
the horizontal stripes those along tp2 and and those along tpi and tp4, which are shown in Fig. 
16. 

In the upper part of Fig. 17, the sequences of ground state solutions for the ola-type structure 
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Fig. 15. The Vp/Vc-dependence of the solution which gives the lowest energy for ^d-type structure, for the case of 
K = 0.25 eV and U = 0.7 eV. 




Fig. 16. Schematic representations of the charge patterns for the different horizontal stripes and vertical stripes in 
the a-type structure: vertical stripes along stack I (a), those along stack II (b), horizontal stripes along and t^z 
(c) and those along tpi and tp4 (d). The thick dotted, thick solid, double dotted, thin dotted bonds represent the 
bonds with transfer integrals tpi, tp2, tp3, tp4, respectively and the thin solid bonds are for tci~3- 



for the case of isotropic intersite Coulomb interactions, Vc = Vp = V, are shown. As V is increased, 
the state which is stable varies from the one with vertical stripes on stack II, to the one with 
another vertical stripes on stack I at V/U = {V/U)cti = 0.27. KF concluded that the latter state 
was the ground state of a-(ET)2l3 for V = under the assumption that the magnetic cell consists 
of four molecules which is identical to the structural unit cell,& but the considerations of larger 
cells in the present calculation indicate that at y = the vertical stripe solution with the other 
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Fig. 17. The sequences of the solution which gives the lowest energy for a-type structures, as a function of V/U 
{U — 0.7 eV). The dotted bonds represent the transfer integral tp4,. The distinction between different stripes among 
the vertical stripes on stack I or II and among the horizontal stripes along tp2 and tp3, and tpi and tp4, are also 
described here. 



charge pattern, i.e. the stripes on stack II, has lower energy. As V is increased further, a phase 
of horizontal stripes along tp2 and tps appears at V/U = {V/U)cr2 = 0.43 and a bipolaronic state 
takes place above V/U = {V/U)cr3 = 0.50. 

In Fig. 18, we show the lp/14-dependence of the ground state solutions for the case of fixed 
value of Vc=0.25 eV. As in the calculations for the 6 and ^^^-type structures, the solution which has 
the lowest energy easily changes from that for the case of isotropic interstite Coulomb interactions, 
when the value of Vp/Vc is deviated from 1. The solution with horizontal stripes along tp2 and tp^ 
has the lowest energy for Vp/Vc < {Vp/Vc)cr = 0.96, whereas the one with vertical stripes on stack 
I is stable for Vp/Vc > iVp/Vc)cr- 

In the work of KF, Hartree-Fock calculations for U = 0.7 on the aMHg-type structure suggest a 
paramagnetic metallic (PM) ground state.& By including Vij some CO phases appear, as can be 
seen in the lower part of Fig. 17 where the ground state solutions for the case of isotropic intersite 
Coulomb interactions, Vp = Vc = V, are shown. The ground state changes from the PM phase, as 
V is increased, to a CO phase with the vertical stripes on stack I at V/U = {V/U)cri = 0.32. As 
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Fig. 18. V^jVc dependence of the solution which gives the lowest energy for the ala-type structure with = 0.25 
eV and V = 0.7 eV. 



V is increased further a phase transition is present at V jlJ = (l^/f/)cr2 = 0.46 to a phase of the 
horizontal stripes along tpi and tp4, and above V jU = {V/U)cr3 = 0.54 the bipolaronic state has the 
lowest energy. As can be seen by comparing the upper and lower parts of Fig. 17, the sequences of 
ground states for the aMHg and als-type structures are qualitatively similar for V/U > {V/U)ct1, 
although the charge pattern for the horizontal stripes stable for {V/U)cr2 < V/U < {V/U)ct3 are 
different to each other. In Fig. 19, the charge density, the magnetic moment on each site, and the 
band gap for the solution with vertical stripes on stack I for the QI3 and oMHg-type structures 
as a function of V/U is plotted which are stable for (y/U)cTi < V/U < {V/U)cr2- The notation of 
'band gap' is the same as in the work of KF,i^ which is the value of the bottom of the 1st band 
minus the top of the 2nd band for the actual case with the unit cell of 4 sites (m = 4 in eq. 3). 
It can be seen that both solutions show qualitatively similar features though the CO state for the 
aMHg-type structure has relatively small values of band gap compared to those for the als-type 
structure. 

3.4 Simplified model for (ET)2X 

As have been seen in § |3.1| - |3.3| the CO states are more stabilized by the intersite Coulomb 
interactions Vij. In contrast, although the presence of relevant values of Vij can be expected in the 
actual K-type compounds which have large dimerization, the dimeric AF state in their insulating 
phase is confirmed experimentally to be stable,& consistent with the MF calculations by KF without 
Vi^j. Here to investigate the competition between this dimeric AF state and the CO states, and to 
get insight into a unified view of the electronic properties of (ET)2^, calculations are performed 
on the model in Fig. 3(b) by varying tpi or tp^. 

Only the case for the isotropic intersite Coulomb interaction V is considered and the values of tc 
and tp are fixed at 0.01 eV and 0.1 eV, respectively. We just take into account a vertical stripe CO 
state described in the inset of Fig. 20. Other stripe-type solutions may have lower energy in some 
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Fig. 19. V/U dependence of the charge density (a), the magnetic moment (b) and the band gap (c) on each site for 
the solutions with vertical stripes on stack 1 for aMHg-type and als-type structures for U = 0.7 eV. 



parameter region but the qualitative features on the competition between the stripe-type CO state 



and other states, whose investigation is our aim, should be made. Note that we know from §3.1 



that this state is stable in the case of tpi = tp^ = tp, identical to the 9-type structure with tp = 0.1 
eV and tc = 0.01 eV, and V/U < 0.49. 

First we fix tp4 = tp, and vary the Coulomb interactions U, V and the degree of dimerization tpi. 
The phase diagram for y = on the plane of U and tpi is shown in Fig. 20(a), which is similar 
to the corresponding phase diagram for the K-type structure obtained by KF^^^ again suggesting 
that this simplified model can be considered as an effective model for the family of (ET)2X. In 
the large U region, the AF insulating state emerges which is the dimeric AF state, and becomes 
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Fig. 20. Phase diagrams of the model structure in Fig. 3 on the plane of U and tpi, in the case of tp4 = tp = 0.1 
eV for 1/ = (a) and for V/U = 0.25 (b). COI, AFI, AFM and PM denote the charge ordered insulating , 
antiferromagnetic insulating, antiferromagnetic metallic and paramagnetic metallic phases, respectively. 



more stabilized as tpi is increased. When V is turned on, there appears the CO phase. The phase 
diagrams on the plane of U and V/U for tpi = 0.30 eV and tpi = 0.15 eV are shown in Fig. 21(a) 
and (b), respectively. It can be seen that the stripe-type CO phase is stabilized when V exceeds 
some critical value, e.g. 0.150 eV for tpi = 0.15 eV and U = 0.7. The bipolaronic state has the 



lowest energy when V is larger than approximately J7/2, similar to the results in §3.1-3.3. If we 
consider the stripe-type CO state with 1 hole on each site along the stripes and the bipolaronic 
state with 2 holes on every 4 sites which will be stabilized in the limit of [/, y ^ oo, the energy 
loss in the Coulomb interactions in the stripe-type CO state and the bipolaronic CO state will be 
2V and U per unit cell of four sites, respectively, so the critical value of y ~ U /2 can be easily 
understood. As can be seen in Fig. 22, where the charge density and the magnitude of spin density 
at each site together with the band gap are displayed for tpi = 0.15 eV, the phase transition from 
the AF insulating phase to the CO insulating phase at V^ri is a second-order transition for U = 0.8 
eV, whereas the one from the metallic phases to the CO insulating phase is a first-order one for 
U = 0.7 eV and 0.6 eV. The stripe-type CO phase which has appeared by the inclusion of V can be 
described on the plane of U and tpi as in Fig. 20(b), where the phase diagram in the case of fixed 
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Fig. 21. Phase diagrams of the model structure in Fig. 3 in the plane of U and V/U, for the case of tpi = 0.30 eV 
(a) and tpi = 0.15 eV (b), where tc and tp4 = tp are fixed at 0.01 eV and 0.1 eV, respectively. The notations for 
COI, AFI, AFM and PM are the same as in Fig. 20, and BI denotes the bipolaronic insulating phase. The dotted 



Unes show the relation 2tpi + (7/2(1 — -v/l + {4tpi/U)^) — V, which is discussed in §4. 4 



ratio of V/U = 0.25 is shown. It can be seen that for tpi < 0.2 eV, where the dimerization is smah, 
the dimeric AF phase is overwhelmed by the stripe-type CO phase, whereas when the dimerization 
is large enough as tpi ^ 0.2 eV this dimeric AF phase survives. 

Next we fix tpi = tp, and vary tp4, so that the band overlap can be controlled, corresponding to 
the case of a- type compounds. The phase diagram for V = on the plane of U and tp4 is shown in 
Fig. 23, where the general features are similar to the Hartree-Fock phase diagram for the a-type 
structure obtained by KF Jli^ as expected. The decrease in the value of tp4 leads to the stabilization 
of the CO state, where the charge pattern and the spin configuration are schematically shown in 
the inset of Fig. 23(b). By the inclusion of V the stripe-type CO phase is stabilized further. Phase 
diagrams on the plane of U and V/U for the case of tp4 = 0.03 eV and 0.09 eV are shown in Fig. 
24(a) and (b), respectively. It can be seen that the stripe-type CO phase is stabilized in a wide 
range of parameters, while the bipolaronic state again has lower energy for V >U /2. In Fig. 25, 
the charge density and the magnitude of spin density at each site together with the band gap are 
shown for the choice of tp^ = 0.09 eV. Here the phase transition as a function of V is of second- 
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Fig. 22. V dependence of the charge density (a), the magnetic moment (b) on each site, and the band gap (c), for 
the case of tpi = 0.15 eV, tp4 = tp = 0.1 eV, tc = 0.01 eV and U = 0.8, 0.7, 0.6 eV. 

order-type when it is from the CO metalhc phase to the CO insulating phase for U = 0.8 eV and 
0.7 eV, and the one from the PM phase to the CO insulating phase is a first-order one for U = 0.6 
eV. The phase diagram on the plane of U and tp4 for the fixed ratio of V/U = 0.25 is shown in Fig. 
23(b), where the CO phase extends compared to the case ofV = 0, especially in the region of large 
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Fig. 23. Phase diagrams of the model structure in Fig. 3 on the plane of U and Jp4, in the case of tpi =tp =0.1 eV 
for V = G {a) and for V/U = 0.25 (b). COI, COM and PM denote the charge ordered insulating, charge ordered 
metallic and paramagnetic metallic phases, respectively. 



§4. Discussions 

In the previous section, it was seen that stripe-type CO states are stabiUzed in general, while 
the actual spatial pattern of charge density sensitively depends on the parameters of models, i.e. 
the anisotropy of the transfer integrals and the values of intersite Coulomb energies. In order to 
find out whether these CO states are realized in the actual compounds, and, if they are, which 
type of stripes are realized, the comparison between our calculation and the properties observed in 
experiments will be pursued in the following, since the value of the intersite Coulomb energies Vij 
as well as its degree of anisotropy Vp/Vc is difficult to determine by theories above. In pursuing 
this, however, we adopt from the studied-* mentioned in §|| the ratio of the neighboring Coulomb 
energies Vc, Vp to the on-site one U to be Vc/U^ Vp/U ~ 0.2 — 0.5 and Vp/Vc to be about 1. 

To reconcile the results of the MF calculations showing ordering of spins to the observed properties 
where magnetically ordered ground states are frequently destroyed by quantum fluctuation which is 
neglected in our approximation, the effects of this quantum fluctuation is incorporated by mapping 
the stripe-type CO states to the 5 = 1/2 Heisenberg Hamiltonians where these spins are interacting 
by the superexchange interactions Jij deduced from the transfer integrals in Fig. 1. Although 
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Fig. 24. Phase diagrams of the model structure in Fig. 3 in the plane of U and V/U, for the case of tp4 — 0.03 eV 
(a) and tp4, = 0.09 eV (b), where tc and tpi — tp eV are fixed at 0.01 eV and 0.1 eV, respectively. The abbreviations 
for COI, COM and PM are same as in Fig. 23, and BI stands for the bipolaronic insulating phase. 



the estimation of the values of Jij along the stripes will be done in the following by using the 
relation Jij ~ "^tij/U with U = 0.7 eV, we note that the definite values should not be taken 
seriously, but only the relative values should, since there is ambiguity on both the value of U and 
the transfer integrals based on the simple semi-emperical extended Hiickel method. However, as a 
general trend, the exchange couplings get reduced from these values of Jij ~ AtfjU in the charge 
ordered state with mtermed.ate charge dispr„p„rt,„nati„„£i thus .t . natural that the exchange 
couplings deduced from experiments are smaller than those naively estimated from this relation. 

4-1 9-phase 

The general features in the phase diagram in Fig. 10 can be understood as follows. In the region 
of small V, the holes tend to form stripes along the direction with smaller transfer energy tc, i.e. 
the vertical stripes, so that they can gain kinetic energy along the bonds with tp. On the other 
hand, the holes are expected to have much localized character in the region of large V, then to gain 
exchange energy along the stripes they choose the stripes with larger exchange energy Jp ~ 4tp/U, 
i.e. the diagonal stripes along tp. 

When tc = 0, the topology of our model is identical to that of the 2D square lattice model with 
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Fig. 25. V dependence of the charge density (a), the magnetic moment (b) on each site, and the band gap (c), for 
tc = 0.01 eV, tpi =tp = 0.1 eV, tp4 = 0.09 eV and U = 0.8, 0.7, 0.6 eV. 



nearest neighbor transfer integral t = tp. The results above for small V that the vertical stripes 
along the smaller transfer integral t = tc are stable, is similar to the results in theoretical studies 
on the Hubbard model on such 2D square lattice, where the stability of stripes along the bonds for 
U/\t\ < 4 and those along the diagonal direction of the bonds for U/\t\ > 4 has been indicated,^!* 
since our calculations were performed for the value of U/\tp\ = 7, which correspond to the latter 
case, and the vertical stripes in our notation are along the diagonal direction of the bonds with the 
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larger transfer integral tp, although the filling of n = f/4 in our calculations is different from the 
above studies for n = 1/2 — x with x ~ 0. 

Prom the MF calculations, the candidates for the ground state of the insulating phase in ac- 
tual 0-type compounds are the CO states with vertical stripes and those with diagonal stripes. 
The bipolaronic state is not likely realized since the ground state in these compounds are mag- 
netic, as disclosed by the magnetic susceptibility measurement .0 When the CO states with the 
vertical stripes and the diagonal stripes are mapped to ID S = 1/2 Heisenberg models, the AF 
couplings between the localized spins can be deduced to be the order of Jc ~ 4t^/C/ = 10 K 
and Jp ~ ^ip/U = 500 K, respectively. The magnetic susceptibilities measured in experiments 
for 6l-(ET)2RbZn(SCN)4 with rapid cooling condition and for 6'-(ET)2CsZn(SCN)4 show Curie-like 
behaviors,0il' suggesting small values of exchange coupling between localized spins, which may in- 
dicate the vertical stripes, also inferred from the reflectivity measurement on 0-(ET)2RbZn(SCN)4 
and 0-(ET)2RbCo(SCN)4 with rapid cooling condition© On the other hand the NMR data on 
0-(ET)2RbZn(SCN)4 with rapid cooling condition suggest that the charge density on molecules be- 
comes continuously distributed by decreasing temperature,© as in the case of the incommensurate 
CDW state. This is in contrast with situation in the ^^-phase discussed in the next subsection and 
apparently inconsistent with our stripe-type CO state suggesting only two kinds of molucules with 
different charge densities. However, since there are two degenerate CO states with vertical stripes 
in the 9-type structure, the dynamical motions of stripes may exist down to low temperature, mak- 
ing difficult to detect the stripe-type CO state by NMR experiments. This charge fluctuation may 
also be the origin for the apparently small couplings between spins, inferred from the magnetic 
susceptibility measurements. The interchain couplings between the spin on these vertical stripes 
show frustration due to the structure (see Fig. 1(b)), which is consistent with the absence of AF 
longe range order down to 2 K.il^ 

Nevertheless, a recent X-ray study on another member of this family, ^-(ET)2CsCo(SCN)4,0^ 
shows a diffusive spot at (0,0,1/2) below 20 K which is proposed to be due to the existence of 
short range CO along the c-axis with period of 2 FT molecules, which is the case of the horizontal 
stripes or the diagonal stripes. The discrepancy between this fact and the conclusion above of the 
existence of CO state with vertical stripes in ^-type compounds based on experimental results on 
members with MM' = RbZn, RbCo and CsZn, may be due to the tendency toward structural 
instability in 6'-(ET)2CsCo(SCN)4, i.e. toward the Od-type structure as in 0-(ET)2RbZn(SCN)4, 
and the fact that X-ray experiments are more sensitive to the lattice distortion than to the charge 
disproportionation. This possibility has also been pointed out by T. Mori et al from thermoelectric 
power measurements.Ei' Another possibility is that the charge pattern of the stripes is different 
between the members of 6'-(ET)2MM'(SCN)4 with MM' = RbZn, RbCo and CsZn, and the salt 
with MM' = CsCo, although their ET layers are isostructural. 
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4-2 9 d- phase 

The stability of the CO state with horizontal stripes in our MF calculations for the ^^-type 
structure compared to the results for the 9-type one where the horizontal stripe solutions have 
rather high energies, may be related to the fact that the PM state in the ^^-type structure, which is 
the ground state for small U and V, shows slight charge disproportionation with pattern of charge 
density similar to that of the horizontal stripe-type CO state described in Fig. 15. 

There are several candidates for the CO state realized in the insulating phase of the actual 
compounds in ^^^-phase based on our calculations. Among them it can be deduced that the vertical 
stripes are not realized in 0rf-(ET)2RbZn(SCN)4 since the temperature dependence of the magnetic 
suscpetibility data for 15 K <T< Tstr = 195 K suggest a larger exchange coupling, J ~ 160K, 
than those inferred from the vertical stripes, Jd ~ 

4t2^/[7 = 50 K and Jc2 ~ ^tl2/U = 100 K. From 
this experimental fact, it is clear that the bipolaronic state with no spins again can be excluded. 
The remaining candidates are the diagonal stripes and the horizontal stripes along tp4 as shown in 
Fig. 15. The CO state with diagonal stripes leads to Heisenberg chains of four kinds of couplings, 
Jpi — JpAi which may have a spin gap at the ground state but inconsistent with the Bonner-Fisher 
behavior in the magnetic susceptibility data for 15 K < T < Tgtr mentioned above. The CO state 
with horizontal stripes along tp4, which we consider to be most likely, is mapped to Heisenberg 
5 = 1/2 chains with a uniform exchange coupling Jp4 ~ ^t^^/U = 450 K, which is consistent with 
this Bonner-Fisher behavior and also with a reflectivity measurement suggesting the existence of 
the horizontal stripes in this compound.ii* If this charge pattern is realized, the spin gap behavior 
seen in experiment is probably due to the spin-Peierls transition along these chains. Note 

that the the lattice distortion expected in the spin-Peierls transition for this case does not produce 
superlattices since the unit cell along the 6-direction initially contains two ET molecules in the 
^d-type structure. The reflectivity data on the ^^-phase of another member (ET)2RbCo(SCN)4 
with slowly cooling condition also suggest the horizontal stripes,® thus it may be concluded that 
the 0rf-type structure in general favors the horizontal stripes. 

Finally we comment on the recent experiment by Miyagawa et aJll' who have performed a careful 
^^C-NMR measurement on 0-(ET)2RbZn(SCN)4 on the condition of slowly cooling, where the 6^- 
type structure is observed below Tgtr- Their conclusion of the existence of two kinds of ETs with 
equal population in the ground state is consistent with our conclusion that the horizontal stripes 
along tp4 are realized in this compound, since in this horizontal stripe solution the charge density 
n{u) on the crystallographically equivalent molecules '1' and '4', and '2' and '3' (see Fig. 1) are 
respectively equal, i.e. n(l) = n(4) and n(2) = n(3), which is not the case in other candidates 
mentioned above showing more than four kinds of ETs with different charge densities {e.g. see 
Fig. 14). Their data showing a broad line and a sharp Pake doublet, ascribed to to the molecules 
with rich and less charge density, respectively, at 15 K < T < T^tv seem to suggest that the charge 
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fluctuation along the stripes are large while in the transverse direction the stripes are 'pinned', 
which may explain why the resistivity shows a sudden increase at Tgtr as lowering the temperature. 
At around the temperature where the broad line turns into a clear Pake doublet (~10 K) the 
magnetic susceptibility also shows the drop. If it is the spin-Peierls transition temperature, this 
result implies that the lattice distortion along the stripes rapidly suppress the charge fluctuation. 

Here we have not discussed the physical origin of the structural phase transition at T = Tstr , since 
the lattice degree of freedom was not included in our calculations. However, it is possible that the 
emergence of the CO state with horizontal stripes and the lattice distortion from 6-type to ^^-type 
structure occur cooperatively, as been discussed in La based high-Tc Cuprates, in the context of 
the relation between the stripe formation and the low temperture tetragonal lattice distortion. 

4.3 a-(ET)2l3 and a-(ET)2MHg(SCN)i 

The difference between results for the 0X3 and aMHg-type structures concerning the charge pat- 
terns in the stable horizontal stripe solution for {V/U)cr2 ^ V/U < (y/U)cr3 can be understood 
when the exchange energies along the stripes, Jpi-pi ~ 'itpi_p^/U, are considered. The two horizon- 
tal stripes are those along tp2 and tp^ and those along tpi and tp4, lead to gains in exchange energies 
of Jp2 ~ 4:tp2/U and Jps ~ Atp^/U, and Jpi ~ Atpi/U and Jp4 ~ Atp/^/U, respectively. In the als- 
type structure, the value of |tp4| is quite small compared to jtpi-sl resulting in Jpi + Jp4 < Jp2 + Jp3, 
thus the calculations showing the stability of the horizontal stripes along tp2 and tp3 with larger 
gain of exchange energies can be explained. On the other hand, in the aMHg-type structure the 
values of Jpi-pA are similar but Jpi -|- Jp4 is slightly larger than Jp2 + Jp3, which is consistent with 
the calculations indicating the horizontal stripes along tpi and tp^ and the calculated energies of 
the two CO states with horizontal stripes are close to each other. 

An explanation we propose here for the observed properties in a-(ET)2l3 is that the stripe-type 
CO state is realized below Tmi, though the charge pattern is different from the one proposed by 
KF which is the vertical stripes on stack 1% The candidates for its ground state based on our MF 
calculations are two kinds of vertical stripe states, a horizontal stripe state and the bipolaronic 
state, as shown in Fig. 17. The bipolaronic state may explain the experimental results since a 
nonmagnetic behavior is seen in the susceptibility measurements, though the required value of 
V/U > {V/U)ct3 = 0.54 in our calculation is rather large. The estimations of exchange coupling 
along the stripes for the stripe-type CO states lead to uniform Jd ~ At'^^/U = 50 K in the case of 
the vertical stripes along stack II, and to alternating Jc2 ~ 41^2/11 = 140 K and Jc3 4t^2/^ = 30 
K in the case of those along stack I. The CO state we propose here to be realized in a-(ET)2l3 is 
the horizontal stripes along tp2 and tp3, with alternating exchange couplings Jp2 ~ ^^^2/^ ~ HOO 
K and Jp3 ~ 'itlr^/U = 220 K, since the observed magnetic susceptibilitylll' below Tmi agree well 
with the calculated susceptibility for the alternating Heisenberg chain with 550 K and 110 K, which 
is respectively half of the the values Jp2 and Jp3, as compared to the one for the Heisenberg chains 
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Fig. 26. Calculated magnetic susceptibilityB^ for the Heisenberg chain models and the measured spin susceptibility 
of Q-(ET)2l3 taken from ref. El 



mapped from the CO state with vertical stripes, as shown in Fig. 26. 
The effect of applying pressure observed in this compound can be inferred from the results 



of calculations on the simplified model in §3.4, as follows. The pressure makes the distances of 
molecules shorter so that the values of transfer integrals increase while it hardly affects the on-site 
Coulomb energy U, so the effective on-site Coulomb interaction i.e. U/W get reduced, where W 
denotes the band width. In Figs. 23 and 24, it can be seen that the PM phase is next to the CO 
phase in the regionof smaller U/W, which is stabilized under high pressure in a-(ET)2l3, as has been 
discussed by KF.& Although the transport measurements by applying pressure and/or magnetic 
field have been intensively done in this compound,0) experiments on the magnetic properties such 
as NMR are crucially needed to understand their electronic properties, since in the CO state the 
charge and the spin degree of freedoms are strongly related to each other. Note that it is difficult 
to decide whether Vij/W increases or decreases by applying pressure to the ET salts since it is 
natural that the values of Vij increase by applying pressure as in the case of the tranfer integrals tij 
and the anisotropy of the ET molecules leads to the difficulty in deducing the the values of intersite 
Coulomb energies Vij between ET molecules, although the values of transfer integrals between ETs 
are carefully calculated by the extended Hiickel method for the cases of different configurations.0^ 
As for Q-(ET)2MHg(SCN)4, the CO metallic phase found in the calculation for the simplified 



model in §3.4 may be the candidate for the ground state supporting the existence of a charge 
disproportionated state proposed by some groups.© However, in §3^ the results on the calculations 
for the actual aMHg-type structure, only insulating states with CO are found to be stable though 
the band gap for V ~ 0.2 has a quite small value of ~ 0.05 eV in the CO state with vertical stripes 
as seen in Fig. 19. So at this point incommensurate SDW state due to the nesting of the Fermi 
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surface, which is not treated in our formulation only considering the states commensurate to the 
lattice, is the candidate for the ground state as discussed in KF.& To understand their electronic 
properties, more studies on this series is needed. 

4-4 Unified view of (ET)2X 

In the work of KF,i)a unified view of (ET)2X was given as shown in Fig. 27(a), where the two 
planes are based on Hartree-Fock calculations considering n and a-type structures on appropriate 



Hubbard models. Based on the calculations in §3.4 showing that the inclusion of intersite Coulomb 



interactions gives rise to a wide region of CO state, their phase diagram can be modified to the one 
as shown in Fig. 27(b), where the two planes are based on the results in Figs. 20(b) and 23(b). The 
location of each type of compound is also indicated in Fig. 27(b): K-phase is in the region of large 
dimerization where the dimer AF state is stabilized, in contrast with 0(i-phase which is located in 
in the region of small dimerization where the CO state is stable. 0-phase is located between al^, 
and aMHg-phases on the axis of band overlap. The main difference in the plane of band overlap 
between Figs. 20(b) and 23(b) is that the aMHg-phase is located in the PM phase far from the 
CO phase in Fig. 27(a) whereas it is close to the boundary between the CO phase and the PM 
phase in Fig. 27(b), so the criticality toward the CO phase may play some crucial role, which is 
beyond the present study. 

The effects of pressure are indicated by arrows in Fig. 27(b), following the discussion of Kanodai^ 
and KF& for the k and a-phases, and that of H. Mori et a for the 6 and ^^^-phases, respectively. 
The directions of arrows for the a and K-phases are opposite from those for the 9 and ^^-phases, 
since in the former compounds it is known that the application of external pressure increases 
the bandwidth, namely reduces the correlation effectj^& while in the latter case it widens (j) and 
results in the decrease of the bandwidth, i) From this point of view, it can be expected that applying 
pressure to the 0-type compounds which is metallic down to low temperature such as ^-(ET)2l3 will 
give birth to an insulating ground state with stripe-type CO state. The effect of uniaxial pressureS^ 
is also an intriguing method since it can tune the anisotropy of the system and the stripe-type CO 
states found in the present study must be highly sensitive to the anisotropy. 

There are discussions that in the case of strongly dimerized structures as in the K-phase the 
effective on-dimer Coulomb energy can be estimated as f/dimer — 2|tpi| + C//2(l — + {^tpi/UY), 
where tpi is the intradimer transfer integral.i'i'il' f/dimer corresponds to the effective on-site Colomb 
energy when one maps models of 1/4-filling with strong dimerization to models with 1/2-filled band, 
e.g as has been studied recently for K-type compoundsil^ for the appropriate 1/2-filled Hubbard 
model. It can be said that C/dimer represent the energy scale in the Mott insulating state. On the 
other hand, as can be seen in Figs. 22 and 25, it appears that the band gap in the CO state increases 
linearly when V is increased, which is similarly seen in a classical CO system Fe304.il'il^ In Fig. 
21(b), the line of C/dimer = ^ is drawn by the dotted line, which shows qualitative feature similar 
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Fig. 27. Schematic phase diagrams for the unifed view of (ET)2X in cases oi V — (a) by KF& and V ^ (b) 
by the present work. The notations AFk and AFq are adopted from ref. ^, which represents the AF phases found 
in their Hartree-Fock calculation on k and a-type structures, respectively. PM and CO denotes the paramagnetic 
metallic and charge ordered phases, respectively. The effect of pressure are shown by arrows. 



to the boundary between the CO state and the dimeric AF state in the case of large dimerization 
tpi = 0.30 eV and in the region of large U. 

§5. Summary and Conclusion 

CO phenomena in 2D organic conductors (ET)2X have been investigated within the Hartree 
MF approximation on relevant extended Hubbard models including both the on-site and the in- 
tersite Coulomb interactions among which the latter are crucial to understand the CO phenomena 
in 1/4-filled systems. By taking into account the explicit anisotropy of the transfer integrals and 
the intersite Coulomb interactions for the 6, 6d and a-type structures, it is found that the inter- 
site Coulomb interactions give rise to stripe- type CO states whose charge patterns are sensitively 
depending on the parameters of these models. 

Based on the results of calculations, implications from the existing experimetal facts lead us that 
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stripes along the bonds with the transfer integral tc are realized in 6-type compounds, those along 
tp4 in Od-type compounds and those along tp2 and tps in a-(ET)2l3. Their magnetic properties are 
well reproduced with these stripe-type CO states, considering the effects of quantum fluctuation by 
mapping the CO states to the S = 1/2 Heisenberg models. A unifed view on the electronic states 
in (ET)2^ is obtained by taking the degree of dimerization and that of the band overlap as key 
ingredients, where several states compete to each other, which are the Mott insulating state, the 
stripe type CO state and the paramagnetic metallic state, realized in the actual compounds. 

Such CO states due to intersite Coulomb interaction are expected to be widely realized in other 
2D non 1/2-filled (non dimerized) organic compounds showing insulating behavior, e.g. a', a" and 
/3"-phases of (ET)2X and also in (ET)mX„ with (m, n) 7^ (2, 1) which do not give 3/4-filled vr-band, 
as well as in non ET compounds.B© The present study, which makes a guideline to understand 
the interrelationship among the experimental facts by taking the CO phenomena as a key factor, 
might help the understandings of their physical properties. 
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